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Abstract 

We demonstrate the existence of a novel breather mode in the self-consistent electron dynamics 
of a semiconductor quantum well. A non-perturbative variational method based on quantum 
hydrodynamics is used to determine the salient features of the electron breather mode. Numerical 
simulations of the time-dependent Wigner-Poisson or Hartree equations are shown to be in excellent 
agreement with our analytical results. For asymmetric quantum wells, a signature of the breather 
mode is observed in the dipole response, which can be detected by standard optical means. 
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Introduction. — The many-electron dynamics in nanoscale semiconductor devices, such 
as quantum wells and quantum dots, has recently attracted a great deal of interest, mainly 
in view of possible applications to the growing field of quantum computing Particular 
attention has been devoted to intersubband transitions, which involve excitation frequencies 
of the order of the terahertz {2]. On this time scale, various collective electronic modes can be 
excited. For instance, the electric dipole response is dominated by a strong resonance at the 
effective plasmon frequency. This resonance (known as the Kohn mode [3() is characterized 
by rigid oscillations of the electron gas, which, for perfectly parabolic confinement, are 
decoupled from the internal degrees of freedom. 

In this paper, we show the existence of a new distinct resonance - a monopole or 
"breather" mode - which corresponds to coherent oscillations of the size of the electron gas 
around a self-consistent equilibrium. Breather modes have been described in many areas of 



physics, such as nuclear matter 



and ultracold atom dynamics 



5| (where they are known as giant monopole resonances) 
7[. In experiments on metallic nanoparticles, monopole 
oscillations of the ionic structure have been observed, which manifest themselves as slow 
modulations of the surface plasmon 8j. However, to the best of our knowledge, previous 
investigations have not addressed the features of the breather mode in the self-consistent 
dynamics of a confined electron gas. Although quantum wells constitute a typical instance 
of such confined systems, the present approach should equally apply to metal nanoparticles 
and carbon-based systems such as fullerenes. 

Model. — Because of the translational symmetry in the transverse plane, the problem re- 

n n 

duces to a one-dimensional (ID) one in the x direction [9|, [10[. To model the electron 
dynamics, we use a self-consistent quantum hydrodynamic model (QHM) that was origi- 



nally derived for quantum plasmas 
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12j and metallic nanostructures 13j. In the QHM, 



the evolution of the electron density n(x,t) and mean velocity u(x,t) is governed by the 
continuity and momentum equations 
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U Tx = -^n^-^^ + ^Yx\-V^^ (2) 
where m* is the effective electron mass, h is the reduced Planck constant, P(x,t) is the 
electron pressure, and V e R = V con f(x) + Vn(x,t) is the effective potential, which is composed 
of a confining and a Hartree term. The Hartree potential obeys the Poisson equation, namely 



Vff = —e 2 n/e, where e is the magnitude of the electron charge and e is the effective dielectric 
permeability of the material. The term proportional to h 2 on the right-hand side of Eq. ([2]) 



12J - or equiv- 



represents the quantum force due to the so-called Bohm potential 

The above QHM can be derived from the self-consistent Hartree equations 
alently from the phase-space Wigner-Poisson equations 15[ - in the limit of long wavelengths 
compared to the Thomas-Fermi screening length. For the sake of simplicity, we shall neglect 
exchange/correlation corrections and assume Boltzmann statistics, which is a reasonable ap- 

nn 

proximation at moderate electron temperatures T |9|, [16J . We also stress that the ID model 
relies on the separation of the transverse and longitudinal directions, which may be broken 
by collisional effects. However, such effects should not be dominant on the fast time scales 
considered here 

The pressure P(x, t) in Eq. (T5]) must be related to the electron density n via an equation 
of state (EOS) in order to close our system of electron fluid equations. We take a polytropic 
relation P = nksT (n/n) 7 , where fee is the Boltzmann constant, 7 = 3 is the ID polytropic 
exponent, and n is a mean electron density. For a homogeneous system (where n = no), 
this EOS correctly reproduces the Bohm-Gross dispersion relation 12j in quantum plasmas. 
For the inhomogeneous electron gas considered here, the choice of n is subtler and will be 
discussed later. 

We assume a parabolic confinement, with V con f = ^u^m^x 2 , where the frequency ujq 
can be related to a fictitious homogeneous positive charge of density uq via the relation 
ojq = (e^o/m^e) 1 / 2 . We then normalize time to cOq 1 ; space to L = (ksT '/m*) 1 / 2 /uq\ 
velocity to L uo; energy to ksT; and the electron number density to n . Quantum effects 
are measured by the dimensionless parameter H = hido/hsT. 

We shall use typical parameters that are representative of semiconductor quantum wells 



10j: the effective electron mass and the effective dielectric permeability are, respectively, 
m* = 0.067m e and e = 13 s$, the equilibrium density is no = 4.7 x 10 22 m -3 , and the filling 
fraction n/no = 0.5. These values yield an effective plasmon energy huo = 8.62 meV, a 
characteristic length Lq = 16.2 nm, a Fermi temperature Tp = 51.8 K, and a typical time 
scale Uq 1 = 76 fs. An electron temperature T = 200 K then corresponds to a value H = 0.5. 

Lagrangian approach. — In order to derive a closed system of differential equations de- 
scribing the breather mode, we fist express the quantum hydrodynamical equations in a 
Lagrangian formalism. We stress that this approach is not based on a perturbative expan- 
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sion, and thus is not restricted to the linear regime. The Lagrangian density corresponding 
to the system of Eqs. ([I])-® reads as (normalized units are used from now on) 
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where the independent fields are taken to be n, 9, and Vh- The velocity field follows from 
the auxiliary function 9(x,t) through u = 89 /dx. The quantity W(n) in Eq. ([3]) originates 
from the pressure, W = f n j^r^f = (3/2) (n/n) 2 . Taking the variational derivatives of the 
action S = J Cdxdt with respect to n, 9, and Vh, we obtain the Eqs. ([I])-©, as well as 
the Poisson equation for Vh- 

The existence of a pertinent variational formalism can be used to derive approximate 
solutions via the time-dependent Rayleigh-Ritz trial-function method 181 ] . For this purpose, 
we assume the electron density to have a Gaussian profile 



n(x, t) 



A 



a 



exp 



(x - df 
2 a 2 



(4) 



where d(t) and a(t) are time- dependent functions that represent the center-of-mass (dipole) 
and the spatial dispersion of the electron gas, respectively. The constant A = N/ \/2ir, is 
related to the total number of electrons in the well, N = j ndx. The above Ansatz is a 
natural one, because for a negligible Hartree energy, V e s reduces to a harmonic oscillator 
potential. 

The other fields to be inserted in the action functional are 9 and Vh- The natural way 
to choose them is by requiring that the continuity and Poisson equations are automatically 
satisfied. The continuity equation is solved with n given by Eq. (00) together with u = 
d + (<r/cr)£, which leads to 9 = (a/2a)£ 2 + d£, where £ = x — d. An irrelevant gauge 
function was discarded in the calculation of 9. The solution of the Poisson equation with a 
Gaussian electron density is 



-A a /2<J — A 
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const, 



(5) 



where Erf is the error function. The integration constant is chosen so that V/f(±a) = 0, 
with 2a being the total size of the system, and letting a — > oo at the end of the calculation. 
As the potential Vh is not bounded, a divergence appears in the Lagrangian density when 



integrated over space. However, the divergent term does not depend on the dynamical 
variables d and a, so that it can be ignored. 

Using the above Ansatz, one obtains the Lagrangian 

1 f „ J d 2 + a 2 d 2 + a 2 
L dx 
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which only depends on two degrees of freedom, namely the dipole d and the variance a. The 
Euler-Lagrange equations corresponding to the Lagrangian L read as 

d + d = (7) 
V2A V3A 2 H 2 
2 3 n 2 a 3 4 a 3 

The quantum-well potential V con { manifests itself in the harmonic forces on the left- 
hand side of both Eqs. flTJ) and (JSJ). As expected, the equations for d and a decouple for 
purely harmonic confinement. Equation (171) describes rigid oscillations of the electron gas 

n n 

at the effective plasmonic frequency, i.e. the Kohn mode [3j, ll9[. Equation (jHJ) describes 
the dynamics of the breather mode, which features coherent oscillations of the width of 
the electron density. The three terms in the right-hand side of Eq. (JHJ) correspond to the 
Coulomb repulsion (Hartree term), the electron pressure, and the quantum Bohm potential, 
respectively. The breather equation (JHJ) can be written as a = —dU/da, where U{a) is 
a pseudo-potential defined as U = o 2 /2 - V2Aa/2 + A 2 / (6 n 2 a 2 ) + H 2 /(8a 2 ). From 
the shape of the pseudo-potential (Fig. H]), it follows that a will always execute nonlinear 
oscillations around a stable fixed point <Jo{A, H), which is a solution of the algebraic equation 
U'(a ) = 0. 

So far, we have not specified the value of the average density n that appears in the EOS, 
P = n 3 /n 2 , written in normalized units. It is natural to assume that n takes some value 
smaller than the peak density at equilibrium Ajo^. The correct way to compute this value 
is to average the square of the density using n itself, i.e. n 2 = (n 2 ) = f n 3 dx/ f ndx = 
^4 2 /(\Z3(Tq). A useful check can be performed by plugging this expression into Eq. (JHJ) 
and neglecting the Hartree potential, which yields the equilibrium variance ero = [1 + (1 + 
H 2 ) 1 / 2 ] 1 / 2 / \J~2. This expression displays the correct low- and high-temperature limits for the 



quantum harmonic oscillator: cr — > 1, for H — > 0; and a ~ a/ H/2, for H 3> 1. 
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FIG. 1: Profiles of the pseudo-potential U(a), for H = 0.5, A = (solid line) and H = 0.5, v4 = 1 
(dashed line). The fixed points are <ro = 1.03 (A = 0) and ao = 1.43 (A = 1). 

With the above prescription for n, the pseudo-potential becomes 
o"q/(2ct 2 ) + H 2 /(8a 2 ). The frequency Q = Q(A,H) of the breather mode, corresponding 
to the oscillations of a, can be obtained by linearizing the equation of motion [Eq. ((SJ)] in 
the vicinity of the stable fixed point of U(a). The dependence of the breather frequency 
with A (i.e., with the electron density) is shown in Fig. [2J For A = (i.e., without the 
Coulomb interaction) the exact frequency is Q = 2uq- For finite A, the breather frequency 
decreases and approaches Q = ujq, for A — > oo. The latter limit can be understood by 
noting that for large A the electron density becomes flatter and flatter, due to the strong 
Coulomb repulsion. Thus, in the limit A — > oo we end up with a uniform electron density 
exactly neutralized by the ion density background. For such a homogeneous system, the 
Bohm-Gross dispersion relation holds, which for long wavelengths yields Q = uq. Indeed, 
if one computes the average density using the prescription used for the EOS, one obtains 
(n) = J n 2 dx/ J ndx = A/(\f2o ) — > 1, for A — > oo (see the inset of Fig. |2]). Thus, as 
expected, for large Coulomb effects the average electron density becomes equal to the ion 
background density. 

Simulations. — In order to check the validity of the above results, we performed numerical 
simulations of the Wig ner-Poisson (WP) system, which is equivalent to the time-dependent 



Hartree equations 15[. In the normalized variables, the Wigner pseudo-probability distri- 



bution f(x,v,t) satisfies the evolution equation 

df df . fdx'dv' rTr i(v-v')x' 

S + 'S- , /SiP f,! "" r "" ,y ' , ' = * (9) 

and is coupled to the Poisson equation. In Eq. (j9|), 8V e s = V e g(x + x' '/2, t) — V e g(x — xf/2, t). 
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FIG. 2: The breather frequency SI as a function of A, for H = 0.5. Solid line: analytical results 
from the Lagrangian method. Squares: The Wigner-Poisson (WP) simulations. The inset shows 
the mean electron density (re) = A/y2ao as a function of A. 

It is important to note that this is a microscopic quantum mean-field model, much more 
general than the hydrodynamic model on which our Lagrangian theory was based. 
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FIG. 3: Solid lines: Electron density profiles at oo^t = 150 from the WP simulations for H = 0.5 
and two values of A. The dashed lines represent Gaussian distributions with the same width and 
the same area as the numerical curves. 

The initial condition used in the simulations is a quantum canonical distribution for the 



harmonic oscillator at finite temperature 20j, where the spatial width <7n has been adjusted 
to the value obtained from the Lagrangian approach to account for the Coulomb repulsion. 
This is very close, but not quite identical, to an exact equilibrium of the WP equations, so 
that the width of the electron density starts to oscillate. We then compute the evolution 
of the dispersion (x 2 ) 1 ^ 2 = (J fx 2 dxdv/ J fdxdv) 1 ^ 2 and its frequency spectrum, which 
generally shows a sharp peak at a dominant frequency. 
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The results of the WP simulations are plotted in Fig. [2] (squares) and agree very well with 
the theoretical curve based on the Lagrangian approach. The agreement slightly deteriorates 
for larger values of A, because the electron density deviates from the Gaussian profile due to 
strong Coulomb repulsion. This is clearly visible in Fig. El where we represent the evolved 
density profiles for two values of A. For A = 1, the profile is still approximately Gaussian, 
whereas for A = 3 an intricate internal structure has developed. Nevertheless, even in this 
case, the error on the frequency is still just over 3%. 

TABLE I: The breather frequency and the equilibrium width <7o for A = 1 and for various values 
of H. 



H o"o (theory) Q (sim.) 

0.00 1.41 1.58 1.60 

0.50 1.43 1.59 1.60 

1.00 1.47 1.60 1.51 

1.50 1.52 1.61 1.57 

2.00 1.59 1.63 1.63 

3.00 1.73 1.66 1.70 



Table [T] shows that the breather frequency depends weakly on the parameter H (and 
hence on the electron temperature). Our theoretical results are in good agreement with the 
simulations, except for H = 1. For this case, the frequency spectrum is particularly broad, 
denoting a significant fragmentation of the resonance. 

Nonparabolic wells. - - For a parabolic potential well, there is no coupling between the 
breather and dipole modes, which may render the experimental detection of the breather 
mode by optical means difficult to realize in practice. 

It can be shown that nonparabolic corrections do not introduce any linear coupling if the 
confining potential is symmetric. A more interesting situation arises for asymmetric wells 



2l|, which we model by adding a small cubic term to the confining potential, V cu t>(x) = 
(K/3)x 3 . The equations of motion are then 

d + d= -K(a 2 + d 2 ), (10) 

V2A y/3A 2 H 2 , . 

a + a = ^— + + —,- 2Kad. 11 

2 3n 2 a 3 4a 3 v ' 



8 



Linearizing (10) and (11) around the stable fixed point (g?o,0"o)> we indeed find a coupling 
between the breather and the dipole, with resonant frequencies fl±. 
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FIG. 4: The frequency spectra of the dipole (solid line) and breather (dashed) modes, for an 
asymmetric well with A = 1, H = 0.5, K = —0.1. 



We model the coupling to the laser field by instantaneously shifting the initial dipole 
of a small quantity, i.e., d(0) = do + d. Figure H] shows a typical spectrum obtained from 
the numerical solution of Eqs. (fTU|) -( lPT|) . which is proportional to the optical absorption 
spectrum commonly measured in the experiments. The sharp peaks correspond to the 
resonant frequencies f2_ = 0.97 and Q + = 1.58, which are rather close to those obtained for 
parabolic confinement. 

The breather mode can thus be triggered using a purely dipolar excitation, and a clear 
signature of the breather frequency can be observed in the optical absorption spectrum. This 
opens the way to optically detecting the breather mode by means of standard pump-probe 
experiments. Finally, the methods used here could be readily extended to 3D nanostructures, 
and may find applications in related areas such as quantum free-electron lasers 22]. 

This work was partially supported by the Alexander von Humboldt Foundation and by 
the Agence Nationale de la Recherche. 
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